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q , Abstract 

We introduce a new spinfoam vertex to be used in models of 4d quantum grav- 
ity based on SU(2) and SO(4) BF theory plus constraints. It can be seen as the 
conventional vertex of SU(2) BF theory, the 15j symbol, in a particular basis con- 
structed using SU(2) coherent states. This basis makes the geometric interpretation 
of the variables transparent: they are the vectors normal to the triangles within each 
tetrahedron. We study the condition under which these states can be considered 
semiclassical, and we show that the semiclassical ones dominate the evaluation of 
quantum correlations. Finally, we describe how the constraints reducing BF to 
\ gravity can be directly written in terms of the new variables, and how the semiclas- 

sicality of the states might improve understanding the correct way to implement 
. ■ the constraints. 

in 
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O ■ 1 Introduction 
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The spinfoam formalism for loop quantum gravity (LQG) [1] is a covariant approach to 
the definition of the dynamics of quantum General Relativity (GR). It provides transition 
amplitudes between spin network states. The most studied example in the literature is the 
Barrett-Crane (BC) model [2]. This model has interesting aspects, such as the inclusion 
of Regge calculus in a precise way, but it can not be considered a complete proposal. In 
particular, recent developments on the semiclassical limit show that it does not give the 
full correct dynamics for the free graviton propagator [3]. In this paper we introduce a 
new model that can be taken as the starting point for the definition of a better behaved 
dynamics. 

Most spinfoam models, including BC, are based on BF theory, a topological theory 
whose relevance for 4d quantum gravity has long been conjectured [4], and has been ex- 
ploited in a number of ways. In constructing a specific model for quantum gravity, there 
is a key difficulty of quantum BF theory that has to be overcome: not all the variables 
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describing a classical geometry turn out to commute. This fact has two important con- 
sequences. The first is that there is in general no classical geometry associated to the 
spin network in the boundary of the spinfoam. This leads immediately to the problem of 
finding semiclassical quantum states that approximate a given classical geometry, in the 
sense in which wave packets or coherent states approximate classical configurations in or- 
dinary quantum theory. This is the problem of defining "coherent states" for LQG, which 
has raised an increasing interest over the last few years [5, 6]. The second consequence 
concerns the definition of the dynamics. This is typically obtained by constraining the BF 
theory, a mechanism well understood classically, but still unsettled at the quantum level. 
The constraints involve non-commuting variables, and the way to properly impose them 
is still an open issue. In particular, it can be argued that the specific procedure leading to 
the BC model imposes them too strongly, a fact which also complicates a proper match 
with the states of the canonical theory (LQG) living on the boundary. For a recent dis- 
cussion of this, see [7]. Notice that these key difficulties are present for both lorentzian 
and euclidean signatures. 

Here we consider a definition of the partition function of SU(2) BF theory on a Regge 
triangulation (for a review see for instance [1]), where the dynamical variables entering the 
sum have a clear semiclassical interpretation: they are the normal vectors n t)T associated 
to triangles t within a tetrahedron r. The classical geometry of this discrete manifold 
I 3X6 clS ; clll gles, volumes, etc.) can be described in a transparent way in these variables, 
provided they satisfy a constraint for each tetrahedron of the triangulation. This con- 
straint is the closure condition, that says that the sum of the four normals associated to 
each tetrahedron must vanish. If this condition is satisfied, this new choice of variables 
positively addresses the issues described above, in the following way. First of all, from 
the boundary point of view, each tetrahedron corresponds to a node of the boundary 
spin network. As we show below, the states associated to the new variables are linear 
superpositions of the conventional ones, with the property of minimizing the uncertainty 
of the non commuting operators: they thus provide a solution to the problem of finding 
coherent states. Upon satisfying the closure condition, the states are coherent and carry 
a given classical geometry. Secondly, as far as the dynamics is concerned, the constraints 
reducing BF theory to GR are functions of the full bivectorial structure of the B field. In 
our model this structure is related to the vectors n t<r in a precise way. We describe how 
this improves the construction of the constraints at the quantum level. We postpone a 
more detailed analysis of the constraints for further work. 

A key point of this approach is the closure condition. As mentioned above, this is 
crucial for the semiclassical interpretation of the states. This is not satisfied by all the 
states entering the partition function. Yet one of the main results of this paper is to show 
that quantum correlations are dominated by the semiclassical states in the large spin 
limit. The key to this mechanism lies in the fact that the coherent states we introduce 
are not normalized, and the configurations maximizing the norm are the ones satisfying 
the closure condition. To prove this, we write the norm as an integral over SU(2), which 
we are able to solve explicitly in particular cases. In the general case, we evaluate the 
integral using the saddle point approximation, and show that the configurations satisfying 
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the closure condition are exponentially dominating. All our calculations are supported 
by numerical simulations, some of which are reported in the Appendix. Remarkably, 
the approximation is accurate even for small spins, and agrees to three digits with the 
numerics at spins of order 100. 

The results we obtain are valid for SU(2) BF theory, but can be straightforwardly 
extended to Spin(4), and the same logic applied to any Lie group, including noncompact 
cases which are of interest for lorentzian signatures. We postpone such developments 
for further work, once the approach considered here proves useful to solve the problems 
described above. 

This paper is organized as follows. In Section 2, we introduce the new vertex ampli- 
tude, and describe how it can be obtained following the conventional spinfoam quantiza- 
tion of BF theory. The new variables are normal vectors associated to the triangles, and 
the new vertex built out of coherent intertwiners. In Section 3 we study the norm of the 
coherent intertwiners, and show that it is exponentially bigger when the closure condi- 
tion holds, thus making semiclassical states dominate the partition function. In Section 
4 we describe the coherent tetrahedron in the conventional basis. In Section 5 we discuss 
how the new vertex can be used as a starting point for models of quantum gravity. We 
summarize our results in Section 6. 

Throughout the paper we use unit h — G — 1. 



2 The new vertex amplitude 

In this Section we describe the vertex amplitude that characterizes our spinfoam model. 
The vertex can be constructed using the conventional procedure for the spinfoam quan- 
tization of SU(2) BF theory. We refer to the literature for details (see for instance [1]), 
and focus here only on the aspect which is relevant to construct the new vertex. 

Recall that there is a vector space H := Inv <^)f =1 Hj i , on which (J^i^i) 2 = 0> 
associated with each edge of the spinfoam. Here the half-integer (spin) j labels the 
irreducible representations (irreps) of SU(2), and F is the number of faces around the 
edge. If the spinfoam is defined on a Regge triangulation, F = 4 for every edge, yet the 
considerations of this paper apply to any value of F. The spinfoam quantization assigns 
to each edge an integral over the group, that is evaluated by inserting in TCq the resolution 
of identity 

!w = \h---3F,il---iF-3)(jl---3F,il---iF-3\- (1) 

The labels % are called intertwiners, and the sums run over all half-integer values allowed 
by the Clebsch-Gordan conditions. Taking into account the combinatorial structure of 
the Regge triangulation (four triangles t in every tetrahedron r, five tetrahedra in every 
4-simplex a), one ends up with the partition function, 

z =ee n^o, (2) 
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where dj — 2j + 1 is the dimension of the irrep j, and the vertex amplitude is A a (j t , i T ) = 
{15j}, a well known object from the recoupling theory of SU(2). 

This partition function endows each 4-simplex with 15 quantum numbers, the ten 
irrep labels j t and the five intertwiner labels i T . Using LQG operators associated with the 
boundary of the 4-simplex, the ten j t can be interpreted as areas of the ten triangles in 
the 4-simplex while the five i T give 3d dihedral angles between triangles (one - out of the 
possible six - for each tetrahedron). On the other hand, the complete characterization of 
a classical geometry on the 3d boundary requires 20 parameters, such as the ten areas plus 
two dihedral angles for each tetrahedron (see next Section for more details). Therefore 
the quantum numbers of the partition function are not enough to characterize a classical 
geometry. 

To overcome these difficulties, we propose to write the same partition function in 
different variables, which will endow each 4-simplex with enough geometric information. 

The key observation is that SU(2) coherent states \ j, h) (here n is a unit vector on the 
two-sphere S 2 - see below and Appendix A for details) provide a (overcomplete) basis 
for the irreps. Given an irrep j, the resolution of the identity in this representation can 
indeed be written as tj = dj f d 2 h \j, h)(j, n\, where d 2 n is the normalized measure on S 2 . 
By group averaging the tensor product of F coherent states, we obtain a vector 

\l,n) : = / dh ®[ =1 h\ji,hi) 

in Ho- Here we denoted j and ft the collections of all jj's and rVs. The set of all these 
vectors when varying the n^s forms an overcomplete basis in H . Using this basis the 
resolution of the identity in TL$ can be written as 
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Ho = / Y[d 2 Ui d k \i,n) {i,n\o = 
J i=i 

p 

= J \\d 2 hi djt J dh J dti h\ji,hi)(ji,hi\ti. (3) 

This formula can be used in the edge integrations instead of (1). The combinatorics is 
the same as before. In particular, notice that (3) assigns to a triangle t a normal n t (r) for 
each tetrahedron sharing t. Within a single 4-simplex, there are two tetrahedra sharing 
a triangle, which we denote u(t) and d(t). Furthermore, (3) assigns two group integrals 
to each tetrahedron, one for each 4-simplex sharing it. Taking also into account the dj 
factors in (3) and using the conventional spinfoam procedure, one gets 

z = e / n rf2 v n d * s n ^> ( 4 ) 

where the new vertex is 

, nt,u(t)\K(t) h d(t) fit ,d(t))- ( 5 ) 
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Let us add a few comments. 

• The new vertex gives a quantization of BF theory in terms of the variables jt and 
n ttT (t)- Together, these represent the full bivector B l (x) discretized on triangles 
belonging to tetrahedra. Taking the B field constant on each tetrahedron, the dis- 
cretization and quantization procedures can be schematically summarized as follows, 

B'Ux) *=5 Bl(r) := Jb^x) dx» A dx» ^ j t n t , T . (6) 

In the quantum theory, the field B l (x) is represented by a vector j t nt, T (t) associated 
to each triangle in a given tetrahedron. The set of all these vectors can be used to 
describe a classical discrete geometry. 

On the other hand, the conventional quantization of BF theory differs in the quan- 
tization step, which reads 

Bl ^ Jl (7) 

where J t are SU(2) generators associated to each triangle. Consequently the vari- 
ables entering the partition function (2) are irrep labels j t and intertwiner labels 
i T . The first variables are related to discretization of the modulus of B, and thus 
to the area of triangles. The intertwiner labels are related to one angle (out of six) 
for each tetrahedron. Therefore in these variables a classical geometry is hidden, 
and this in turn makes it hard to implement the dynamics, as we discuss below in 
Section 5. 

• This vertex can be used directly for SU(2) BF theory, and straightforwardly gener- 
alized to the Spin(4) case, exploiting the homomorphism Spin(4) = SU(2) x SU(2). 
In particular, (6) becomes 

Bit(x) ^ B { j (t) := Jb%(x) dx» A dx» ^ j+ n+ r © j~ n~ T . (8) 

• There is a simple relation between this vertex and the one given in (2) which is the 
one usually used. Since all we did is simply a change of basis, (5) can be written as 
a linear combination of {15j}'s as we will see in more details in Section 4. 

• The BC vertex can be obtained from (5): performing the d 2 hi integrations with the 
condition h ttU{t) = h tA{t) (namely inserting terms S(nt, u (t) ~ nt,d(t)) in (4)), we indeed 
obtain the BC vertex in its integral representation [8], 

a bc = [ n ^ n t x {jt \Kt t) h m ). (9) 
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To understand the geometry of this new vertex amplitude, it is crucial to study its 
asymptotic behaviour in the large spin limit. 1 In the rest of this paper, we begin the 
analysis of this limit, by studying in details the behavior of the states \j,n)o entering (5). 
The full asymptotics of the vertex will be reported elsewhere. Nonetheless, it will become 
clear from the rest of the paper that the large spin limit of this integral is dominated by 
values of the group elements such that the factors of the integrand in (5) are close to one, 
namely such that 

\jt, n tM t)) = K(t) h d{t) \jt, n tAt) ). (10) 

This has a compelling geometric interpretation. Recall that nt, u (t) and n tt d(t) are the 
normals to the same triangle as as seen from the two tetrahedra sharing it. Then (10) 
means that they are changed into one another by a gravitational holonomy. 

This is in contrast with the BC model, which fixes n ttU M = nt,d(t), thus not allowing 
any gravitational parallel transport between tetrahedra. This is another way of seeing 
the well known problem that the BC model has not enough degrees of freedom. We will 
come back to this point in Section 5, where we discuss how this vertex can be taken as 
the starting point for quantum gravity models. 

3 Coherent intertwiners 

In this Section we focus on the building blocks of the new vertex, namely the states \ j,n)o 
associated to the tetrahedra entering (3). Before studying the details of the mathematical 
structure of these states, let us discuss their physical meaning. In the canonical picture, a 
tetrahedron is dual to a 4-valent node in a spin network. More in general, when the edge 
of the spinfoam is on the boundary, there is a one-to-one correspondence between the 
number of faces F around it, and the valence V of the node of the boundary spin network. 
Given a discrete atom of space dual to the node with V faces, its classical geometry can be 
entirely determined using the V areas and 2(V — 3) angles between them. Alternatively, 
one can use the (non-unit) normal vectors Hi associated with the faces, constrained to 
close, namely to satisfy J2i ™i = 0- 

On the other hand, the conventional basis used in (1) gives quantum numbers for V 
areas but only V — 3 angles. This is immediate from the fact that one associates SU(2) 

— * — » — » 

generators J« with each face [9], and only V — 3 of the possible scalar products Jj • Jk 
commute among each other. Thus half the classical angles are missing. To solve this 
problem, we argue that the states \ j,n)o carry enough information to describe a classical 
geometry for the discrete atom of space dual to the node. In particular, we interpret 
the vectors j% precisely with the meaning of normal vectors n« to the triangle, and we 
read the geometry off them as described above. This requires the vectors to satisfy the 

1 The reasons why the large spin limit is relevant to discuss the dynamics of spinfoams is largely 
discussed in the literature (see for instance [1]). The main reason is that this limit appears to be related 
to Regge calculus, a discrete version to classical GR. 
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following constraint: 

v 

n = J2^^ = °- ( n ) 

i=i 

This closure condition will be studied in detail in the rest of this Section. If (11) holds, 
all classical geometric observables can be parametrized as 0({jihi}). 

This is the physical picture we consider. To describe the mathematical details, let us 
recall basic properties of the SU(2) coherent states (for details see the Appendix A). A 
coherent state is obtained via the group action on the highest weight state, 

\j,n)=g(n)\j,j), (12) 

where g(n) is a group element rotating the north pole z = (0, 0, 1) to the unit vector n (and 
such that its rotation axis is orthogonal to z). These coherent states are semiclassical in the 
sense that they localize the direction n of the angular momentum. We have (j, n\J\j, n) = 
j h and relative uncertainty A ~ 1/ y/j vanishing in the large spin limit (see Appendix 
B). As \ j,j) has direction z with minimal uncertainty, \j, h) has direction n with minimal 
uncertainty. 

This localization property is preserved by the tensor product of irreps <8>i\ji,fii): in 
the large spin limit we have for instance (Jj • Jk) — jijkhi ■ hk for any % and k. More in 
general any operator 0{{ Ji }) on the tensor product of coherent states satisfies 

(6({Ji}))*o(tii*i}) (is) 

with vanishing relative uncertainty. 

So far, the ?Vs are completely free parameters. If the closure condition (11) holds, 
(13) tells us that the states admit a semiclassical interpretation. Consequently, states 
in <S)^ =1 'Hj i satisfying (11) can be considered semiclassical. Yet they cannot be used in 
the partition function: the tensor product of coherent states must be projected on the 
invariant subspace Ho, in order to implement gauge invariance. This can be achieved by 
group averaging, as anticipated in the previous Section, and the resulting states written 
as a linear combination of the conventional basis of intertwiners are 

H 3 \hQ) := dh®Y =1 h\j i ,n i )= Q 1 ...j v _ 3 (^, n*) \j ± . . . j v ,ii ■ ■ ■ ivs)- (14) 

h...i V - 3 

We call these states coherent intertwiners. The explicit form of the Cj 1 ...j v _ 3 coefficients 
is reported in the Appendix A, but it will not be relevant for the rest of the paper. In 
Section 4 below we will study their asymptotics in the 4-valent case. 

Crucially, (13) holds on the states \j,n)o (see Appendix B), thus group averaging does 
not spoil the semiclassical interpretation of these states in the large spin limit. But what 
is more important, group averaging correlates the ?Vs, making the state sensitive to (11). 
Indeed, notice that a state in Ti.o satisfies ( Yli Ji) = not necessarily the closure 
condition. The main result of this paper is to show that the states satisfying (11) have 
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an exponentially bigger norm in the large spin limit. Therefore every state (14) enters 
the partition function (4), but quantum correlations will be dominated by semiclassical 
states in the large spin limit. We can then say that projecting onto Ho does not force 
the ^-simplex to close identically, but it does imply that closed simplices dominate the 
dynamics. 



3.1 The closure condition 



The key for the mechanism is the following: the coherent interwiner states (14) are not 
normalized, and the norm is a function of the hi, 



f{hi) := o{j,n\j,n) = / dh (ji, hi\h® v \j h hi) ^ 1. 



(15) 



Consequently, the relative weight of these states will depend on the configuration of the 
hi. As we show below, the norm is exponentially smaller when the closure condition (11) 
is not satisfied. 

To give a simple picture of the way this happens, let us first consider the trivial case 
of a bivalent node. In this case, the norm (15) can be straightforwardly evaluated using 
the explicit values of the Clebsch-Gordan coefficients. Using \ji,hi) = g%\ji,ji), we can 
write 

(j^m,^) = D^(g-^D^(h)D^( 9i ). (16) 

The integration over h in (15) gives the orthogonality condition for the Clebsch-Gordan 
coefficients, 



J dhD$\h)DW(h) 



(17) 



and we get f{h\,fi2) 



This matrix element can be evaluated 



starting from the fundamental representation. As we show in the Appendix B, we have 



D^ j M 1 92)\ 2 =\(M9i 1 92\-)[ il = ( 



1 - - {hi + h 2 



> 2 ) 



2jl 



(18) 



Defining J = we conclude that 



f{hi,h 2 ) 



§31 32 



d 



1 - 



31 



AT2 

IP 



2jl 



(19) 



As N 2 < J 2 , we see that in the large spin limit the norm vanishes exponentially, unless the 
closure condition is satisfied, in which case it vanishes with an inverse power law. Con- 
sequently, the norm of the states satisfying the closure condition becomes exponentially 
bigger than the others. 

Notice that in the bivalent case, the closure condition means that the two vectors 
are antiparallel. The situation is analogous for the trivalent case: closed configurations 
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will correspond to three coplanar vectors. It is only for valence four or higher that the 
geometry of the problem becomes more interesting, and one has closed configurations 
spanning 3d space. However, the configurations where the vectors are all aligned are 
still present. Since such configurations fail to produce 3d structures, we refer to them as 
degenerate. 

For nodes with higher valence we expect a result analogous to (19). However the 
analysis is far more intricate, the reason being that now when we use (16) we obtain an 
integral J dh D^ ni (h) of V representation matrices, and in general we have no simple 
formula for the (generalized) Clebsch-Gordan coefficients (see Appendix A) appearing 
from this integration. It is then simpler to study the norm explicitly using the fact that 
(see Appendix B) 

/ 1 1 \ 2 i 

{j,j\g(n)~ 1 hg(n)\j,j)= +\g(n)~ 1 hg(n)\-, +)J = (00S7 + ism-fu ■ h) 2j , (20) 

where we parametrized SU(2) group elements as 

h(j,u) = cos 7 1 + i sin7-u • a, 7 £ [0, 7r], ueS 2 . (21) 

It is convenient to introduce the vector p = sin7w. In these variables, the Haar measure 
splits into two integrals over the unit three-ball B : \p\ < 1, corresponding to the northern 
(say B + ) and southern (say B ) "hemispheres" of SU(2) ~ S 3 , 




We then write (15) as 

/<".) = h e/ b -A rife- + >p- < 22 > 

v =± JB v V 1 P i=l 

where p v = rj^/l — p 2 in B n . It is a triple integral depending on 3V parameters. Notice 
that the measure tends to locate the integral on the boundary \p\ = 1 of the integration 
domain, whereas the integrand tends to locate it on the origin \p\ = 0. In the large 
spin limit, the latter behaviour will of course dominate. Since we are interested in the 
asymptotic behavior of this norm, it is enough to study its saddle point approximation in 
the large spin limit. In the Appendix C we discuss the exact evaluation of the norm. 



3.2 Saddle point analysis 

We want to study the behavior of (22) in the large spin limit. To this end, we scale the 
spin labels by an overall constant A, i.e. ji 1 — > A A;^ and write (22) as 

/(*) = 2^£ Sb ~7T^ eXm ' S(?):=£2fe In (p,, + (23) 



r,=± 
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Notice that the real part of S is negative. When the variable A is taken to infinity, the 
integral can be approximated computing the saddle point expansion and evaluating the 
Gaussian integral on the maxima of S. 2 It is immediate to see that the norm of S is 
maximized by p = 0, corresponding to h = 1, but the phase requires more attention. The 
saddle points of S(p) on B v satisfy 

*f P^ + ip-rii) 2 Pv i Pit + ip-ni) 2 

A solution of this set of six real equations can be found taking the scalar product of the 
real part with p, giving 

( 1 ~pJTW 7 w) =0 - 

Since < p 2 + {p ■ hi) 2 < 1 and h L > for all i's, this implies that p 2 + {p ■ hi) 2 = 1 for 
all i, namely p 2 = {p • hi) 2 . 

This leaves us with two cases: 

• Either there exists (at least) two non-collinear vectors, hi ^ then the latter 
equality forces p = and the imaginary part of (24) simply reads J2 i kih; L = 0. 

• Or all the hi are equal to the same unit vector h up to a sign, then a little work on 
the imaginary part of (24) leads us to the same constraint £V fcj hi = but leaves 
no constraint on the vector p. 

In both cases, we see that we have derived the closure condition (11) from the saddle point 
analysis in the large spin limit A — > oo. Numerical investigations show that there are no 
other solutions to (24), and that there are no saddle points if the closure condition is not 
satisfied. Which as we show below means that as A increases the norm (15) is exponentially 
smaller for non-closed configurations, and thus correlations on spin networks in the large 
spin regime will be dominated by those intertwiners satisfying the closure condition and 
thus describing a classical geometry. 

Let us start by the non-degenerate generic situation with non-collinear ?Vs. Then 
S(Q) = and the Hessian matrix of the second derivatives at the two (one for each three- 
ball) fixed points p = is a sum of projectors P l ab = (5 a b — (hi) a (hi)bj orthogonal to the 
unit vectors hi, 

1 " E ^ ( 25 ) 

p=0 i 



H a b = — ~ 



2 dp a dp b 



Notice that this is independent of 77, thus both saddle points give equal contributions. 
Next, the property J2ij hi kj ( n « ' n j) 2 < kj guarantees that the eigenvalues of H ab 



2 Let us point out that the measure divergence at \p\ = 1 does not contribute at all to the leading order 
behaviour of the norm f{hi) unlike the case of the asymptotics of the 6j and lOj symbols [10]. Indeed, 

the divergence in 1 is still integrable and does not lead to a divergent integral. 

yi-p 2 
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are all positive. Therefore the Gaussian integral on the real line converges and we can 
straightforwardly compute the saddle point approximation, 

ftnA = -L / d 3 p e~ XH ° bPaPb = \ , (26) 

We can compute det if explicitly in terms of the labels hi and normals nf. 

K . 1 - 

detH = —J2 k i kj \hi A hj\ 2 - - ^ A;, /c fc \hi ■ (hj A n fc ) | 2 , 

with X = ki. When the closure constraint is satisfied, i.e. £V kihi = 0, the quadratic 
term is a function of the (squared) area of the internal parallelograms while the cubic 
term relates to the (squared) volume of the tetrahedron. From this point of view, the 
determinant det if measures the shape of the tetrahedron at fixed triangle area fcj. 

The result (26) has been numerically confirmed for various closed configurations and 
different values of V. As an example, a plot is reported in Appendix D, which shows 
explicitly how accurate the approximation is, even for very small spins. In particular, we 
have one digit of accuracy from the beginning, and three digits at spins of order 100. We 
conclude that the norm of coherent intertwiners satisfying the closure condition scales as 
\~ 3 / 2 in the large spin limit. 

3.3 Degenerate configurations 

Let us now look at the case where all the h^s are collinear. Without loss of generality, 
we can take all the n^'s aligned with the z axis, namely hi = rjiZ, rji = ±. Let us 
also define ^ = J2i\ Vi =±ji- ^ ne can easily convince himself that 3 /faz) = unless 
Wi^i — 1 and j + = j~ = -j A. Then, (20) reads (cos 7 + rji sm^u z ) 2ji . We take the 
conventional parametrization S 2 3 u = (cos a sin /3, sin a sin /3, cos /3), with normalized 
measure d 2 u = -^dadf3 sin/5. We see that the a angle drops out, and the integral is 
simply 

1 r 

- d(3d^ sin (3 sin 2 7 (1 - sin 2 /3 sin 2 7 ) yA . (27) 

71" Jo 

This integral can be exactly evaluated as shown in the Appendix C, and the result is 

(vx + iy 1 . 

We conclude that 

r vt+i if n,^ = i, j + = r = |a, 

f(\Tk2) = \ (28) 
[ otherwise. 

Notice that the degenerate configurations scale as A -1 , whereas the closed ones scale 
as A~ 3 / 2 . For A ^ 00, the degenerate configurations dominate as y/\. However, in the 



3 It suffices to decompose the state as in (14) and recall the conditions for a tensor product of states 
to admit singlet irreps. 
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partition function (4) they have zero measure, thus we do not expect them to affect 
the dynamics in a relevant way. Furthermore these degenerate configurations are also 
present in BC, and they do not enter the computation of physically relevant quantities, 
as conjectured in [11] and shown in [12]. 

3.4 Non-closed configurations 

— * 

For configurations such that iV ^ 0, numerical simulations show that there are no saddle 
points, but the integrand in (22) is still maximized at p — 0. For notational consistency, 
let us also rescale N i— > XN. The expansion of S in (23) around p = reads 

S(p) = 2tN-p-H abPaPb + o(p 3 ). (29) 

If we approximate the integral using this expansion for the exponent, we see the presence 
of a phase term coming from the linear term in (29). This phase term dumps the integral 
exponentially. In fact, a simple calculation gives 

ffa) = J_ / d ^ e i2XN-p-XH abPaPb = 1 *_ = / -N-H-^Y (3Q) 

This result is again supported by numerical simulations for various non-closed configu- 
rations and different values of V (see Appendix D). As N ■ H~ 1 N > 0, the norm of 
non-closed configurations is exponentially smaller than the closed ones. 

We do not have a simple explicit expression for the inverse Hessian H -1 , but we can 
easily write it as a power series. Notice that H in (25) reads as K(t — E) with K = ^ fc, 
and £ a & = ^Ziiki/K) h\. As the norm of E a fe is smaller than 1, the inverse series for H^ 1 
converges (except in the degenerate configurations where some of the hi are collinear): 

h- 1 = ^[i + e+e 2 + ...] = ^[n-E|<^+E^^-%)n>;+...]. 



This allows us to express the exponent of (30) in terms of the matrix = v ^ 3 (fii-hj), 
N-H^N = j^y^^j kj (hj ■ hj) + kj kj k k (hj ■ hj) (hj ■ h k ) + . . . 

hi i,i,k 

= ^'T;i(; ■ a- ■ a" ■ ...), l ^ r k~. 

Notice that GV,- is a V x V Gram matrix of scalar products of unit vectors. For V greater 
than 3, the vectors hi are not linearly independent and detG^ = 0. Nevertheless the 
series G + G 2 + G 3 + . . . can be re-summed as G(l — G)^ 1 . Thus it appears that the 
speed of convergence of the intertwiner norm in this non-closed case is directly related to 
the eigenvalues of the Gram matrix G. More precisely, estimating its smallest and largest 
eigenvalues (both between and 1) would allow to bound the convergence speed of the 
norm as A goes to infinity. 
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4 Four-valent case: the coherent tetrahedron 



Above we have constructed coherent intertwiners for nodes of generic valence. The 4- 
valent case, whose dual geometric picture is a tetrahedron, is of particular interest as 
it enters the construction of vertex amplitudes for spinfoam models. In this Section, 
we study how the coherent tetrahedron can be decomposed in the conventional basis of 

— * — * 

virtual links, which we fix by choosing to add J\ and J 2 first. 

The coefficients entering (14) can be studied with the same techniques used for the 
norms. Introducing the shorthand notation \i) = \ ji,hi), in the 4-valent case we have 

\(ji---j4,Ji2\j_,n)o\ 2 = j dhj dgd jl2 x {jl2 \g) (12|%|12> (34|/i|34>. (31) 

As it can be verified directly, dj 12 x^ 12 Kq) projects on the intertwiner state \ji . . .j^jvi)- 
To study the asymptotics, it is convenient to introduce an auxiliary unit vector h, writing 
the character in the basis of coherent states, 

X Ul2 \g) = d jl2 [ d 2 n (j 12 ,n\g\j 12 ,n). (32) 



As in the previous Section, we parametrize the group elements as h = p v + ip ■ a, g = 
q r f + iq- <J. We have 

gh = p v q v , -p-q + i(p v q + q v >p-pAq^) - a. (33) 
Using this and (20), and rescaling again ji \— > A ki, we rewrite (31) as 

\(n...H,Ju\j,n)^ = ^-J2 f ^=£=f -f£=[d*ne^ (34) 

with 

/ \ 2fc 

e S(M) = {q^+iq-nf^W^q^ - p ■ q + i ( Pri q + q^ p - p A q) ■ 

i=i 

4 



From the analysis of the previous Sections, we expect the asymptotics of these coefficients 
to be dominated by g and h close to the identity. Expanding S around p = q = and 
denoting r — (p,q), we have 

S(r) = 2iN-r-r ■ Hr + o(r 3 ), (36) 
where N is the following 6-dimensional vector, 
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The Hessian matrix has the following structure, 



I' //' 



with H the same as in (25), and 



F ab = k x Pl b + k 2 P 2 ah -irjrf t abc (\ (ni) c + fc 2 (n 2 ) c ) , (39) 
H' ah = hP^ + hP^ + hiPa,. (40) 

The antisymmetric part of F comes from the pAq term in (33), so from the non-abelian 
nature of SU(2). Notice that the only imaginary term in H is proportional to 7777', thus 
when we perform the sums in the Gaussian approximation, we are simply going to get 
(four times) the real part of a single Gaussian integration. 
The asymptotics are given by 

where det H can expressed in term of 3 x 3 determinants as detif det(G — F T H~ 1 F). 
This integral can be studied numerically. For fixed ji, n,, it represents the probability of 
the eigenstate ju as a function of the n^s. For closed configurations, we expect this to 
be a Gaussian peaked on the semiclassical value computed from the n^s. Let us consider 
for simplicity the equilateral case. In this case, we expect ju to be peaked around j such 

that j(j + 1) = (( Ji + J 2 f) = 2j'o(f jo + 1), which gives J = ^2j (|jo + 1) + J - \- As 
we show in Fig.l, for large spin (41) is indeed approximated by the Gaussian 

PO'12) = N(j ) exp (- (J12 ~ J)2 ) , (42) 



where N(jo) is the normalization. Confronted with the numerics, we fix a = jo/2. This 
semiclassical property is manifestly preserved by changing the pairing, namely the choice 
of basis. For any choice of jik the probability (41) is peaked 01a j ik (j ik + l) ~ (jihi+jkn k ) 2 
with vanishing relative uncertainty. 4 

Analogous results hold for arbitrary closed configurations. This shows in a concrete 
way in which sense \j,n)o represents a semiclassical state for a quantum tetrahedron, and 
more in general for a l^-simplex dual to a node of valence V. 

At this point, it is useful to compare our construction of the semiclassical tetrahedron 
with the one introduced in [6]. Both states have the property that for any i and k, (Ji ■ J k ) 



4 An alternative definition of coherent states, considered in [13], is to require the minimization of the 
uncertainty A J, L -J k A J v J t > ±|([AJ; • J k ,£\J t ■ J k ])\ = \\(e abc J? J h k Jf)\ - O(f). The states constructed 
here do not stricly minimize this, even if the do satisfy the non trivial condition AJi ■ J k AJi ■ Ji ~ 0(j 3 ) 
which guarantees the vanishing relative uncertainties. 
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Figure 1: Plots of (41) for the equilateral configuration ji = jo Vi. The dots represent the exact 
numerical evaluation, whereas the line is the Gaussian (42). The left panel shows the case jo = 10, 
whereas the right panel shows the case jo = 100. In the small spin case, the Gaussian approximation is 
already capturing the right behavior, and it becomes very accurate in the large spin case. 

is peaked around a given semiclassical value with vanishing relative uncertainty in the large 
spin limit. The uncertainties are slightly different, as can be seen by direct comparison 
in the equilateral case. From equation (41) of [6] (with the appropriate redefinition of jo) 
we read a = 4jo/3\/3 which is roughly one and a half times the spread of (42). The main 
advantage of the states introduced here lies in the fact that they allow the construction 
of a wider framework. The construction in [6] is only valid for a tetrahedron, and tailored 
to studying the large spin limit only. The framework considered here, on the other hand, 
applies to any V^-simplex, and it is well-defined for any spin. In particular, the key feature 
of the coherent intertwiners considered here is that they provided an overcomplete basis 
in Hq, which makes them useful to study the dynamics, as we described in Section 2. 



5 Towards the quantum gravity amplitude 

Let us finally discuss how the coherent state technology presented here can be used to 
define the dynamics of quantum GR, starting from the spinfoam model (4). Spinfoam 
quantization of GR usually relies on reformulating GR as a constrained BF theory with 
an action of the following type [14, 15, 16], 

S G n(B^,u„) = J T[BAF(u)+C(B^). (43) 

is a connection valued on a given Lie algebra (for euclidean signature, su(2) or spin(4)), 
and B^ v is a bivector field (or two-form) with values in the same Lie algebra. The term 
C(B) includes polynomial constraints, reducing topological BF to GR. 5 Typically, it gives 

5 The logic of this reduction is the following. The initial BF action describes a topological theory with 
no local degree of freedom, with the field B a Lagrange multiplier enforcing that the connection is flat, 
F(uo) = 0. The term C{B) then constrains the Lagrange multiplier, thus enlarging the phase space. It 
breaks the topological invariance (translation symmetry of B) and introduces non-trivial local degrees of 
freedom. 
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the set of second class constraints expressing B in terms of the tetrad field (or vierbein) 
and leading to GR in the first order formalism (for a detailed canonical analysis of (43), 
see [17]). 

The BF theory can be quantized as described in Section 2, discretizing the spacetime 
manifold with a Regge triangulation (or more generally a cellular decomposition), and 
then evaluating the partition function (2), where the variables are the representations j t 
and intertwiners i T . The natural extension of this procedure to quantize GR with action 
(43) would be to discretize the constraints C(B) and include them in the computation of 
the discretized path integral (see e.g. [18, 19]). Nevertheless, on the grounds of geometric 
quantization (see e.g. [2]), one usually shortcuts this computation and directly attempts 
to translate the discretized constraints C(B) as constraints on the variables j t and i T by 
assuming that the variables B are represented in the spinfoam as the generators of the 
considered Lie group, as in (7) for SU(2). This procedure is usually referred to as imposing 
the constraints at the quantum level, and it leads to a unique choice of intertwiner and to 
the Barrett-Crane model [2, 20]. However, the quantum constraints obtained in this way 
are fairly strong, and a number of authors have raised various doubts on this procedure, 
such as: the rewriting of the constraints in purely algebraic terms (representations and 
intertwiners) partially hides their geometric meaning, and thus the interpretation of the 
model; the unique intertwiner does not seem to be compatible with the data on the 
boundary spin network as given in LQG; more crucially, the quantum constraints do not 
commute with each other and generate (by commutator) higher order constraints which 
do not seem to have any classical equivalent. 

A plausible explanation of these difficulties is that the BC model imposes the con- 
straints too strongly, and thus does not have enough degrees of freedom to describe 
quantum GR. Then a way out would be to loosen the implementation of the constraints, 
by requiring them to hold only on average (vanishing expectation value) with minimal 
uncertainty. This seems the natural procedure when dealing with non-commuting con- 
straints. 6 It has been suggested for instance in [18], where both possibilities are discussed, 
(i) to impose the constraints exactly by inserting the 5-distribution 5(C(B)) in the path 
integral, or (ii) to regularize the ^-distribution by inserting its Gaussian approximation 
exp(— C(B) 2 / £ 2 ) , with £ a free parameter. Notice that the Gaussian insertion amounts to 
impose the constraints approximatively. More precisely, since the B field represents the 
geometry and large values of B corresponds to large (length) scales, using the Gaussian 
means imposing the constraints at the semiclassical level while allowing slight deviations 
at smaller scales. 

In the same spirit, here we argue that the partition function (4) for the quantum BF 
theory in terms of coherent states offers a natural way to impose the constraints on aver- 
age. The key is that the partition function (4) provides a quantization of BF theory where 
the B field is represented not through generators of the group, but as representation and 
intertwiner labels, such as jt and n t (see (6) or (8)). This representation allows to write 
directly the constraints C(B) as a sum of local contributions C a (j t n tT ), in terms of the 

6 A possible strategy to impose strongly second class constraints is discussed in [21]. 
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dynamical variables and with a clear geometric interpretation. In particular, the semiclas- 
sical dynamics becomes transparent. As we understand from the analysis reported in the 
previous Sections, the vertex is dominated in the large spin limit by semiclassical states 
satisfying the closure condition for each tetrahedron and the relation (10) between adja- 
cent tetrahedra in the same 4-simplex. On these states the variables j t and n t give classical 
values with an (almost minimal) uncertainty decreasing as the jVs increase. Therefore, 
imposing C a (j t n tjT ) = amounts to impose the classical constraints only on average with 
a small uncertainty vanishing in the large spin limit, namely semiclassically. If we use the 
ji's and nj's to construct a Regge geometry, we expect that the role of the constraints 
is to generate deficit angles when we glue together various 4-simplices, thus allowing the 
geometry to be curved. We postpone the precise analysis of the implementation of the 
constraints for further work. 

6 Conclusions 

The standard intertwiner basis which leads to BF theory with the vertex amplitude given 
by the { 1 5 j } symbol does not appear to be the most suitable one to study the semiclas- 
sical geometry of BF theory. Furthermore, it makes it hard to understand the quantum 
structure of the constraints reducing BF to GR, and possibly hides the correct way to 
implement them. 

Here we considered a basis constructed out of SU(2) coherent states. We defined non- 
normalized coherent intertwiners, and studied their norm as a function of the geometric 
configuration described. For each configuration, the norm is an integral over SU(2) that 
can be solved exactly as described in the Appendix. Yet a saddle point evaluation of the 
leading order of this integral in the large spin limit proves a very accurate approximation 
even for small spins, and shows very neatly that the norm is exponentially maximized 
by the states admitting a semiclassical interpretation, namely the ones whose quantum 
numbers can be interpreted as vectors ji hi describing the classical discrete geometry of 
a V-simplex. Thanks to this crucial result, the semiclassical states will dominate the 
evaluation of quantum correlations. 

Using these coherent intertwiners we rewrote the BF partition function with a new 
vertex amplitude, given in (5), where the discrete B t {r) variables are interpreted in terms 
of the vectors jt ht >T , thus retaining the original vectorial nature of the B field. We expect 
this reformulation of the BF spinfoam amplitudes to improve the geometric interpretation 
of the theory and in particular our understanding of what should be the proper way to 
implement the constraints reducing it to GR. 

We hope that the developments presented here will contribute to the promising recent 
advances in understanding the low energy limit of LQG, such as the ones on the graviton 
propagator [3, 11, 12, 22], on the emergence of effective actions for matter [23], or on 
applications to cosmology [24] and black hole physics [25] . 
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A Coherent States for SU(2): a brief review 

SU(2) coherent states minimize the (SU(2) invariant) uncertainty A = |(J / — {J) 2 \ 
in the direction of the angular momentum [27]. A simple calculation shows that on a 
basis state \j,m) the uncertainty A(J,m) = j(j + 1) — m 2 is minimal when m = j. 
The maximal weight vectors are thus coherent states for arbitrary choice of spin 
j and angular momentum axis J z . Starting from the highest weight, an infinite set of 
coherent states on the sphere SU(2)/U(1) ~ S 2 are constructed through the group action, 
\j,h) = g(n)\j,j), where n is a unit vector defining a direction on the sphere S 2 and 
g(n) a SU(2) group element rotating the direction z = (0, 0, 1) into the direction n. 

— * 

Explicitly, taking n = (sin 9 cos 0, sin 9 sin 0, cos 9), we choose g(n) = exp{i9 rh ■ J} where 
rh = (sin0, — cos0, 0) is a unit vector orthogonal to both the directions z and n. Just 
as \j,j) has direction z with minimal uncertainty, \j,n) has direction n with minimal 
uncertainty, as can be verified explicitly using the formula reported in the next Appendix. 
A coherent state can be decomposed in the usual basis as 

3 

= ^2 Om(n)\j,m), (44) 

m=-j 

where 

a (h)-< ! m C ~ tan 6 e -» 

These states are normalized but not orthogonal, the scalar product between two of 
them being 

(3M3M) = ( i± 4p) J <»«»^\ A(z,n 1 ,h 2 ) = -'-In (l±||) . 

Here ^ = ^e~' t( ^. Notice that A is the area of the geodesic triangle on the sphere S 2 by 
the north pole direction z and the two unit vectors h\ and hi- 

Consequently they provide an overcomplete basis for the irreps j, and the resolution of 
the identity can be written as tj = dj f d 2 h \j,n)(j,n\, with d 2 h the normalized measure 
on S 2 . 

The explicit coefficients entering (14) can be found using (44) to decompose \j,n) 
into the conventional basis of 7Y , 

v 

n 3 \j_,n ) = II a ^ J2 W I-? 1 ■■■3v.ii--- iv-z) , (45) 

mi. ..my i=l ii...iv—s 
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where we have introduced the (generalized) Clebsch-Gordan coefficients Cmi'.'.'.my (see 
[28]). These are defined 7 by the integration of irrep matrices, 

j d 9 X[D^ ni {g)= C;:r;X : (46) 

i=l ii...ivs 

Using the recoupling theory, thes generalized coefficients and can always be decomposed 
into sums of products of conventional (3-valent) Clebsch-Gordan coefficients. 
From (45) we immediately read the coefficients entering (14). 



B Useful formulas 



We report here the results used in Section 3. The calculations can be done using the 
explicit expression (44), however it is usually easier to exploit the fact that 

(j,n\J a \j,n) = (j,j\J' a \j,j) 

where J' a = gin)" 1 J a g{h) is the rotated generator. 

We compute the following averages of the SU(2) generators on the coherent states: 

j i 

(j, n\J a \j, h) =jn a , (j, h\J 2 \j, n) = - + j(j - -) n a . (47) 
From this it is easy to check that 

O — »Q — * — * 

A = (j,n\J \j,n) - (j,n\J\j,n){j,n\J\j,n) = j. 
Further calculations give 

j 1 

(j,n\J a J b \j,n) = - (5 ab + ie abc n c ) + j(j - -)n a n b , 



(121(7! + J 2 ) 2 | 12) = ( 3l h 1 +j 2 n 2 ) 2 + 3l +j 2 , 
(here notice that the ji + j 2 term comes from the uncertainty A 2 computed above), and 
(12|(J 1 + J 2 ) 4 |12) = ((j 1 n 1 +j 2 n 2 ) 2 + 3l +j 2 ) 2 + 



+2jlj2 



(ji + 32 - ^) (l - (ni • n 2 ) 2 ) + (1 - 3 ni • n 2 ) 



Therefore the uncertainty ((J x + J 2 ) A ) — {{J\ + J 2 ) 2 ) 2 depends explicitly on the value of 
the angle cos#i 2 = hi ■ h 2 . For large ji,j 2 , it is for aligned vectors and maximal for hi 
and h 2 orthogonal. 

Proceeding as above, one can show (13) for any observable O(Ji). 

Let us now report other results used in the main body of the paper. 



7 Here we loosely use the term Cebsch-Gordan coefficients, to refer to invariant tensors. These can 
differ in phase and normalization from other definitions found in the literature. 
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The matrix element entering the bivalent norm (18) can be computed directly using 
the conventional parametrization 

g{hi) = g(9i,rhi) = cos 6 { 1 + i sin Q { rhi ■ a, (48) 

we have 

9i 1 g2\— )| 2 = cos 2 6\ + cos 2 9 2 — 2 cos^i cos 02 (cos 6*i cos 6*2 + sin^i sin 9 2 rh\ • m 2 ) = 

= i -^ = 1 _i ( . i+ft2)2 m 

In Section 3 we also made use of the adjoint action of the group on itself. Using 
(21) and (48), we have 

g{n)~ x hg(n) = cos7 + i siii'ju' ■ a, u' = (?(n) _1 u. 

By definition g(h) rotates the north pole vector (0,0,1) to the direction n, thus 
(■u') 3 = (g(h)~ 1 u) 3 — u • fi, as can be explicitly checked writing the components of 
the rotated vector u'\ 

Q 

u' = cos 9u — sin 9 u Am + 2 sin 2 - (u ■ m) m, 

(■u') 3 = cos 9u 3 — sin 9 (u A m) 3 = cos 9u 3 + sin 9 (ill cos + m 2 sin 0) = -u • h. 
From this we immediately derive (20). 

A crucial result concerns the extension of (13) to expectation values in Ti Q . To prove 
this, one can compute 

(j, n\h J \j, n) — j[ cos n + i sin7-u + sin 7 u A fi\ (cos 7 + i sin7M • n) 2 - 7 ' -1 . 

As we see, the presence of the group element h complicates the expectation values 
of the SU(2) generators (compare the above with (47)). However, we know from the 
analysis in Section 3 that the integration over h peaks the group elements on the 
identity 7 = in the large spin limit, thus 



,n\hJ\j,n) ~jh, 

and all the results reported above extend naturally to Hq in the large spin limit. 
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C Evaluating the norm using spherical integrals 

In Section 3 we wrote the norm (22) using p to parametrize SU(2). This choice was 
convenient to study the saddle point approximation. On the other hand, to evaluate the 
norm exactly, it is more useful the standard parametrization (21) in terms of a rotation 
angle 7 and its rotation axis u, 

f(hi) = / dh(ji,hi\h m \ji,hi) = - sin 2 7 / d 2 u TT(cos7 + i sin 7 u ■ hi) 2j \ 

J n Jo Js 2 f = 7 

We can expand all the terms and deal separately with the integrals over 7 and the spherical 
integrals over u. Using the binomial expansion and J = J2iJii the ^7 integrals are of the 
following type, 

with k < J. 

The next step is to compute the spherical integrals. One can show that: 

2 1 
d 2 u JJ(m -hi) = -ni • h 2 , 



3 



4 



/ 

/ d 2 u TT(m • hi) = [(ni • h 2 ) (n 3 • n 4 ) + (ni • n 3 ) (n 2 • h 4 ) + (n x • n 4 ) (h 2 ■ h 3 )] . 

J i=1 * x 

This can be generalized to even polynomials of arbitrary high order in u. A generic term 
/ nf^C^'^j) wm have a prefactor (2iV+l)!! and a sum over the (N — l)\\ possible pairings 
of the vectors hi with each other. 

Using these formulas one can compute exactly the coherent intertwiner norm f(hi) as 
a finite sum of polynomials of the hi. 

The situation is particularly simple in the degenerate case, when all the hi are aligned. 
In this case, discussed in Section 3.3, the spherical integrals are trivial, and we have to 
evaluate simply 

f(N,r)iz) = ^-Jdpd-f sin (3 sin 2 7(1 - sin 2 p sin 2 ~f) N = (51) 

= ^B" 1 )* ( N k ) I ^( sin ^) 2fc+1 / ^(sin 7 ) 2fc+2 = £ ( N 

where in the last step we used 

J ^(sm/3) 2fc+1 = v^^|, / rf7(sin7) 2fe+2 = v^ r(A; + 1 



r(* + §) 
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This sum in (51) can be straightforwardly evaluated with the substitution t = s + 1, to 
obtain 



f(N n-z)-— T ( - 1} " 1(iV + 1)! - --±-((1 - - l) = - J- (52) 
This result was used to evaluate (27). 

D Numerics 

All the calculations of this paper have been supported by numerical simulations, per- 
formed with Mathematica™. In this Appendix we report some examples to illustrate the 
numerical support to the saddle point analysis of Sections 3.2 and 3.4. 

The first thing we show is the absence of saddle points when the closure is not 
satisfied. Consider the integrand in (22), with the parametrization (21) where u = 
(cos a sin (3, sin a sin (3, cos (3). To give a clear picture of the integrand, we fix an arbi- 
trary value of a and show a 3d plot in only (3 and 7. We choose the simplest case, the 
equilateral tetrahedron (V = 4) with ji = 100 Vi. This is shown on the left panel of Fig.2. 
It is clearly visible that 7 = is a saddle point (the degeneracy in (3 is an artifact of 
the polar coordinates for S 3 ), and that the integrand is quickly suppressed away from it. 
On the right panel, we have again V = 4, but this time we picked a random non closed 
configuration, with the four different spins still of order 100, but not equal. First of all, 
one can notice that even if 7 = still maximizes the integrand, it is not anymore a saddle 
point. Second, notice that the integrand has now also negative values, which contribute 
to dump the integral in this non closed case. 




Figure 2: The real part of the integrand for V — 4 as a function of 7 £ [0, 0.2] (symmetric at tt/2) and 
j3 £ [0, 7r], for fixed a. Left panel: the equilateral case, at ji = 100 for all i. Right panel: a generic open 
configuration, with different spins but all of order 100. 

Next, we consider the Gaussian evaluation of the norm. In Fig. 3 we compare the 
exact numerical evaluations of (22) with the analytic calculations. On the left panel, we 
consider a closed configuration, for simplicity given by the equilateral tetrahedron, namely 
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V = 4 and ji = j Vz. The dots represent the numerical evaluations for different values 
of j, whereas the line is the analytic result (26). As one can see, the agreement is very 
good also at small spins. By direct reading of the numerics, one sees that the analytic 
approximation matches to one digit the exact result already at j ~ 1, and that at j ~ 100 
the matching is to three digits. The situation is very similar for arbitrary configurations. 
On the right panel we consider a non closed configuration with ji = j The dots 
are again the numeric evaluations, and the line the analytic result (30). Notice that the 
agreement is still very good even at small spins. As expected, the norm is exponentially 
smaller than in the closed configuration. In the open case, the integral converges poorly, 
for this reason we have less points in the plot. 




1 2 5 10 20 50 20 40 60 80 



Figure 3: Bilogaritmic plots. Left panel: the dots are the numerical evaluation of the exact norm (15) 
for the equilateral tetrahedron (Ji = j for all i), for different values of the spin j. The line is the analytic 
calculation of the leading order (26). Right panel: an open case with ji — j for all i but the normals not 
closing. The line is the analytic result (30). 

The results reported above show how accurate the saddle point approximation is. 
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